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Abstract 


I 

□ 


In  this  paper ,  we  will  give  an  h-vetsion  finite  dement  method  for  a  two  dimensimial  nonlinear 
dasto-plasticity  problem.  A  family  of  admissible  onistitutive  laws  based  on  the  so-called  gauge, 
function  method  is  introduced  first ,  and  then  a  high  order  h-version  semi-discretisation  scheme  is  _ 
presented .  The  existence  and  uniqueness  of  the  solution  for  the  semi-discrete  problem  are  guaranteed 
by  using  some  special  properties  of  the  constitutive  law ,  and  finally  we  will  show  that  as  the  maximum 
element  sise  h  0  ,  the  solution  of  the  soni-discrete  problem  will  converge  to  the  solution  of  the 
continuous  problem .  The  high  order  h-version  discretisation  scheme  introduced  here  is  unusual .  If 
the  partition  of  the  spatial  space  only  has  rectangles  or  paralldograms  involved ,  then  there  would  not : 
he  any  limit  on  the  element  degree .  However ,  if  the  partition  of  the  spatial  space  has  some  triangular 
elements ,  then  only  certain  combinations  of  finite  dement  spaces  for  displacement  and  stress  functions  - 
can  be  used .  The  discretisation  scheme  also  provides  useful  idea  for  applications  of  hp-version  or  high 
order  h-version  finite  element  methods  for  two  dimensional  problems  where  the  elasto-plastic  body  is  ^ 
not  a  polygon ,  such  as  a  disk  or  an  annulus .  i  - 


Dist 


Introduction 


Several  papers  about  finite  element  methods  with  some  theoretical  convergence  analysi 


Avail  and/or 
Special 


results  for 


two  dimensional  dasto-plasticity  problems  have  been  published  in  the  en^neering  literature ,  such  as  the 


work  of  Bonnetier  [3] ,  Johnson  [10]  [11] ,  and  Miyewhi  [14] .  However ,  in  all  these  approaches ,  the  shape 
of  the  domain  is  assumed  to  be  polygonal  and  is  partitioned  into  several  triangles .  The  displacement 


functions  are  approximated  by  piecewise  linear  functions ,  and  the  stress  functions  are  approximated 
by  piecewise  constant  functions .  This  choice  of  finite  dement  spaces  for  the  displacement  and  stress 
functions  makes  it  easy  to  use  the  constitutive  law ,  because  both  the  strain  and  stress  functions  will  be 


piecewise  constant ,  and  each  element  can  therefore  be  defined  as  in  dther  its  dastic  or  plastic  state . 
However ,  for  smooth  solutions  or  over  sub-re^ons  where  the  solution  is  smooth,  higher  order  polynomial 
spaces  are  usually  preferred .  Moreover ,  if  the  domain  under  consideration  has  a  special  shape ,  such  as  a 
rectangular  domain  or  an  £  shaped  domdn ,  then  we  can  just  use  rectangular  elements  in  the  partition . 
For  rectangular  elements ,  we  have  to  apply  piecewise  bilinear  approximation  for  the  displacement  function 
in  order  to  obtain  a  conforming  finite  dement  method .  In  this  case ,  the  corresponding  strain  functions  are 
piecewise  linear  instead  of  piecewise  constant .  Hence  the  stress  function  can  no  longer  be  approximated 
by  piecewise  constant  functions ,  and  we  can  no  longer  define  a  whole  element  as  elastic  or  plastic .  As 
we  will  sec ,  the  design  of  a  scheme  to  choose  an  appropriate  constitutive  law  on  high  order  elements  is 
an  important  and  subtle  problem . 

’Partially  supported  by  NSF  grant  DMS-91-20877  &  ONR  grant. 

^Partially  supported  by  NSF  grant  DMS-91-20877 . 
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A  family  of  admissible  constitutive  laws  based  on  the  so-cidled  gauge  function  method  was  first 
introduced  by  Bonnetier  [3] .  An  h-version  method  with  piecewise  linear  displacement  function  for  prob¬ 
lems  based  on  a  constitutive  law  of  this  kind  has  been  presented  in  his  work .  This  paper  presents  a  high 
order  h-version  finite  element  method  for  elasto-plasticity  problem  based  on  an  admissible  constitutive 
law .  An  unusual  method  to  define  the  constitutive  law  over  each  element  is  proposed .  The  convergence 
analysis  shows  that  the  limit  function  is  guaranteed  to  satisfy  the  constitutive  law  over  the  whole  domain . 
Two  types  of  elements ,  namely  rectangular  and  triangular  elements ,  are  discussed .  In  addition ,  some 
possible  generalizations  of  our  method  are  also  mentioned . 


2  A  Family  of  Admissible  Constitutive  Laws  for  Two  Dimensional 
Problems 


The  family  of  constitutive  laws  that  will  be  introduced  below  is  based  on  the  following  two  basic 
assumptions : 

1 )  Existence  of  a  conv^  yield  surface . 

2 )  The  normality  condition :  the  plastic  increment  is  proportional  to  the  outward  normal  to 
the  yield  surface  during  plastic  flow . 

There  are  two  main  reasons  why  we  choose  this  family  of  constitutive  laws .  First ,  they  are  actually 
a  generalization  of  some  of  the  most  commonly  used  engineering  formulations .  Second  ,  Bonnetier  has 
shown  that  the  continuous  problem  based  on  this  fatmily  of  constitutive  laws  will  be  well  posed . 

For  two  dimensional  problems,  we  will  describe  a  yield  surface  by  the  stress  tensor  ass  (oji ,  022, 012) 
and  a  set  of  hardening  parameters  ( sometimes  also  called  internal  parameters  ) : 

a  =  ( oi ,  02 , . . . ,  am  C  W  C  II"* 
where  W  is  a  convex  set  in  R™  . 


The  elastic  set  is  assumed  to  be  convex ,  and  hence  we  can  think  of  a  function  F*  ( o ,  a  )  which  is  a 
convex  function  of  a  and  a  as  the  gauge  function  of  this  convex  set .  More  precisely  we  assume  that 
there  exists  a  function  F  :  R^  x  W  • — ►  R  such  that 


A1 ).  F{a,  a)  is  a  convex ,  and  piecewise  amalytic  function  of  o ,  a .  (2.1) 


A2). 

/13). 


F(0,0)  =  0. 


There  are  constants 
uniformly  on  the  set 


dF  dF 

7 ,  r  >  0  such  that  7  <  I  I  ’  I  ^  I  <  ^ 

{(o,a)|F{<T,o)  =  Zo)  for  some 


Zo>0. 


where  ^ 
oa 


f  dF  dF  dF  Y 
\  daw  ’  da22  *  dai2  ) 


and 


dF 

da 


(  dF  dF  dF  Y 
\da^'  da^'""  '  da^) 


(2.2) 

(2.3) 
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The  corresponding  constitutive  law  can  then  be  derived  by  using  the  two  postulates  mentioned 
earlier.  Let  (<r ,a)  lie  inside  the  yield  surface,  i.e.,  /’(<r,a)  <  Zq  .  Then  the  material  is  assumed  to 
be  dastic .  Therefore  the  stress  and  strain  rates  must  satisfy  Hooke’s  law ,  and  there  is  no  change  in  the 
hardoiing  parameters.  So 

f  d  =  Dc 

(2.4) 


I 


d  =  Dc 
d  =  0 


whoe  d  =  ~  ,  e  s  ^  etc ,  and  D  is  the  usual  elasticity  tensor .  In  plane  stress  problems ,  for 
ot  in 

instance ,  the  elasticity  tensor  has  the  form : 

\  V  0 
V  \  0 

0  0  — 


D  = 


l_„a 


where  E  and  v  are  the  Young’s  modulus  and  the  Poisson’s  ratio  of  the  material . 

If  the  point  ( <r ,  o )  lies  on  the  yield  surface  when  <  s  to  t  i-e.  F  (a ,  a)  (to)  ==  Zq  ,  and  will  move 
towards  the  inside  of  the  yield  surface  afterwards,  then  we  stiU  say  that  the  material  is  elastic  at  t  s  to  . 
So  the  stress ,  strain  and  hardening  parameter  rates  will  still  satisfy  the  equation  (2.4) . 

If  the  point  ( ^ ,  a )  lies  on  the  yield  surface  when  t  =  to  ,  and  will  remain  on  the  yield  surface 
for  to  <  t  <  to  +  j  ,  i.e.  f  ( o  (t) ,  a  (t) )  s  Zo «  V  t  €  ( <0 « fo  +  ^ ) ,  then  we  say  that  the  material  is 
plastic  at  t  s  to  and  we  will  split  the  strain  rate  i  into  elastic  and  plastic  parts 


€  =  +  ^, 

and  assume  that  the  elastic  part  is  still  related  to  the  stress  rate  by  the  Hooke’s  law 

d  = 

Then  by  the  normality  condition  3  A  >  0  such  that 


(2.5) 


(2.6) 


r  -p  1 

■  ar  • 

—  \ 

da 

—  d 

•  ^ 

dF 

.  da  - 

(2.7) 


During  the  plastic  flow ,  ( <7(t)  ,o(0 )  still  remain  on  the  yield  surface,  therefore  differentiation  of 
the  equation  F{<r{t)  ,a{t))  =  Zq  with  respect  to  time  i  yields 


dF'^  .  dF"^  .  „ 


(2.8) 
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Using  (2.6)  ,  (2.7)  and  (2.8),  we  can  eliminate  to  obtain: 


OSS  D _ III _ gg  gg  _ _  i 

OF  ’ 

da  da  da  da  • 

< 

dF^. 

da  da 

Hence ,  the  constitutive  equations  based  on  the  gauge  function  f'  (  o ,  a  )  read 


(  a  ^  Di 
I  d  =  0 


for  ia,a)  €  £ 


a  -  fn-  ■ 

^  ’  daF'^dc.F  +  d^F'^DdaF  ^ 


.  _  d^F^a  ft  p 


for  (a, a)  €  V 


where 


=  |(®,o 


)  F(a,a)<Zo  or  F(<r,a)  s  Zo  and 


dF^  \ 


f  ffpT 

’  ss  j  (<^»<»)  F(a  ,a)s:  Zo  and  a>0 


(2.10) 


(2.11) 


(2.12) 


(2.13) 


The  three  assumptions  Al)'vA3)  on  the  gauge  function  F(a,a)  are  essential  for  the  proof 
of  the  existence  and  uniqueness  of  the  solution .  Actually  from  the  assumption  A3 ) ,  we  can  get  the 
following : 

Proposition  2.1  £ei  a  s=  .4(a,a,c)  ,  to/iere  A  is  the  constitutive  operator  defined  by  (2. JO) 
and  (2.11) .  Assume  that  the  gauge  function  F{a,a)  satisfies  the  assumption  AS) .  Then  there  exist 
constants  0  <  7*  <  F’  independent  of  the  values  of  a  jO  and  i  such  that  for  any  (a  ,a)  €JRP  xlf 
and  any  c  £  II"*  , 


7'|c|^  <  ^(a,o,c)^f  <  r'|c|^ 


(2.14) 


Proof  :  Any  vector  c  =  ( cn ,  €22  ,€12)^*  can  be  decomposed  into  a  sum  of  two  vectors ,  one  parallel 
to  d^F  ,  and  the  other  orthogonal  to  daF  ,  with  respect  to  the  scalar  product  {a  ,  b)  =  (Da  ,b)  = 
a'^Db,  i.e. 


e  =  +  €  " 


(2.15) 
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c  *  =  <  d^F,d»F  >-*  <  deF,i>  d,F 

=  {d^F'^Dd^F)-'id„F^Di)d,F 


(2.16) 


iDi^,i')  =  0. 


(2.17) 


Since  D  is  positive  definite  and  independent  of  the  values  of  o ,  a  and  c  ,  it  is  sufficient  to  show 
that  there  are  constants  0  <  7'  <  F*  such  that 


(2.18) 


Now  suppose  the  material  is  elastic ,  then  the  constitutive  operator  A  is  defined  by  (2.10),  and  we  get 


^(<y,o,0^€  *  (^<.<)  =  («.«)• 


(2.19) 


Otherwise ,  if  the  material  is  plastic ,  the  constitutive  operator  A  will  be  defined  by  (2.11)  and  we  can 
let 

n  =  [  daF’^daF  +  daF'^D  ]  '^ 

and  get 

Aia,a,i)^i  =  [{D  --  nD8„Fd„F'^D)i,e] 

=:  {€-  ti{d„F'^Di)d„F,i) 


From  the  definition  of  c  *  ,  wc  have 

A{c,a,ifi  =  (i^,<'-)  +  (l-/.(ftF’'£>»,/’))(i',£') 

By  virtue  of  the  assumption  A3) ,  the  term  udaF^daF  is  bounded  from  below  and  above: 


(2.20) 


n(i  +  l|x?||) 


<  ft(daF^d^F)  <  1 


* 

which  .shows  that  (2.14)  holds  for  y'  =  ^ 


lli^ll) 


and  r’  =  1  . 


Definition  2.1  The  family  of  the  constitutive  laws  (2.10)^  hosed  on  the  gauge  function  of  the 
yield  surface,  F{a  ,a)  which  satisfies  (2.1)'^  (2,S)  will  be  called  admissible  constitutive  laws  or  briefly 
admissible  laws . 


As  we  mentioned  earlier ,  the  admissible  constitutive  laws  are  actually  a  generalization  of  some  of 
the  most  commonly  adopted  engineering  formulations.  For  instance,  the  kinematic  hardening  law  for 
two  dimensional  problem  ( also  called  Ziegler's  rule  ) ,  can  be  formulated  as  follows : 
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Lei  0  s  (/Jj]  ,022%0i2)^  the  hardening  parameters,  and  t;  >  0  be  the  constant  coefficient  of 
hardening  rate .  The  Von  Mises  surface  is  chosen  to  be  the  yield  surface ,  and  it  has  the  form : 

(2.21) 

with  2o  >  0  a  constant  and 

/(ff)  =  (ffii  +  -  oiiOM  +  3o?a  . 

The  constitutive  equations  based  on  Ziegler’s  rule  then  take  the  form : 

{a  —  Di  if  /(o  -  0)  <  Zo 

(2.22) 

0  =  0  or  f(o  -  0)  =  Zq  and  <  0 


Now  let  R  be  the  matrix 

'  4/3  2/3  0 

R  =  2/3  4/3  0 

.0  0  1/3 

and  v/R  be  one  of  its  square  roots .  We  can  then  define  another  set  of  hardening  parameters  a  = 
( y/rjfy/K)~^  0  ,  and  a  gauge  function  of  the  form : 

F(<T,a)  =  f{<r  -  y^y/Ra).  (2.25) 

Obviously  the  gauge  function  F  ( <7 ,  a )  defined  by  (2.25)  is  a  convex  and  analytic  function  of  cr 
and  a  ,  and  satisfies  the  equation  F(0,0)  =  0.  A  simple  computation  also  shows  that  on  the  set 
{  (0,0-)  6  R®  X  R®  I  F(<T,a)  =  Zo  }  ,  the  derivative  d^F  satisfies  the  inequalities 

I  <  \d,F\  <  Vi. 

Since  R  is  positive  definite,  doF  =  -  y/rfy/Kd^F  is  also  bounded  uniformly  on  this  set .  Hence  the 
gauge  function  F  satisfies  all  three  of  our  assumptions .  It  is  not  difficult  to  check  that  in  (2.10)  (2.13) , 

if  we  use  the  gauge  function  F  defined  as  in  (2.25) ,  we  will  then  get  a  set  of  constitutive  equations  which  is 
equivalent  to  the  equations  (2.22)  and  (2.23)  after  the  variable  substitution  a  =  ( y/lj  y/K  )~^  0  .  Hence 
we  can  consider  the  kinematic  hardening  law  as  a  particular  case  of  our  general  admissible  constitutive 
laws . 
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Similarly  the  two  dimensional  pure  isotropic  hardening  law  can  be  cast  into  our  admissible  form  by 
using  a  gauge  function  of  the  form : 


F(a,a)  s=  /(<y)  -  h(a), 


whoe  Or  €  2/  C  R  is  the  isotropic  hardening  parameter  and  h  (  a )  is  a  concave  monotonically  increasing 
function  of  the  variable  a  with  h  (0)  =  0  . 

3  Two  Dimensional  Quasi-static  Elasto-Plastic  Problems  with  Con¬ 
stitutive  Law  Based  on  A  Gauge  Function  Method 

A  typical  two  dimensional  quasi-static 
eiasto-plastic  problem  with  the  constitutive 
law  based  on  the  gauge  function  method  can 
be  formulated  as  follows : 

Let  SI  ,  a  bounded  domain  in 
(cf.  Figure  3.1),  be  the  reference  config¬ 
uration  of  an  eiasto-plastic  body.  We  de¬ 
note  by  U  =  {U\  ,Ui)  the  displacement 
field  of  the  body  where  IJ\  (* ,  y)  i  Ui  (x ,  y) 
are  the  displacements  at  the  point  (z ,  y) 
in  the  x  and  y  directions  respectively .  Let 
a  s  (oii  be  the  stress  tensor, 

where  <T|  |  and  022  arc  the  normal  stress  in  the  x  and  y  directions ,  and  012  is  the  shear  stress .  We 
will  use  e  =  ( cn  ,(22  ><12  to  denote  the  strun  tensor ,  where  the  normal  strains  cn  and  €22  and 
the  shear  strain  e\2  are  defined  in  the  normal  engineering  formula; 


Figure  3.1  A  bounded  domain  in  R^ 


ai/i  dV2  dU^dV2  .... 

‘>'  =  17  ’  '"  =  17  •  =  ‘  ’ 

Finally  I  =  (0,r]  will  be  the  time  interval  of  interest ,  and  /  =  (/i,/2)^  is  the  body  traction  acting 
on  the  eiasto-plastic  body  SI .  Then  a  general  two  dimensional  quasi-static  eiasto-plastic  problem  with 
the  constitutive  law  based  on  the  gauge  function  method  is  as  follows . 

To  find  the  displacement  functions  U{x,y,t)  and  the  stress  functions  a(i,y,t)  ,  such  that 
a{x  ,y,t)  satisfies  the  following  equilibrium  equations : 


a.e.  on  fl  x  / 


(3.2) 
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with  the  boundary  condition : 


Ui{x,y,t)  s  Uiix  yjf.t)  =  0  for  (z,y)  €  dfi, 

and  the  stress  function  o  ( z ,  y ,  t )  defined  by  the  following  constitutive  equations : 

=  Di  (x,y,t) 

tt(*,y,0  =  0 

f  d(z,y,0  =  (2?  -  !)')€(*. y.t) 

a(x,y,i)=  (z , y , 0 


(3.3) 


! 


for  (o,a)  €  2^ 


.(3.4) 


(3.5) 


where  D '  is  defined  as 


nV.  «  -  {Dd.FdeF'^D){x,y,i) 

^  ^  (d«f’7daF  + )(z,y,o  ■ 


(3.6) 


(3.7) 


The  stress ,  strain  and  hardening  parameter  functions  satisfy  the  following  initial  conditions : 

=Q  ;  f/,(z,y,0)  =  0  i,jf=l,2  , 

o,(z,y,0)=0  i=l,2,...,m  . 

Of  course  here  we  have  to  require  that  the  body  traction  functions  satisfy  the  compatibility  conditions : 

/<(*»y»0)=0.  »*1,2  .  (3.8) 

Finally  E  ,  V  ase  the  elastic  and  the  plastic  sets  defined  by : 


2:  =  I  (<T,o)  €  xR." 
F  =  I  (<7,a)  €  R^  xR" 


if 


if 


F(ff,a)<Zo  or 

F(a,o)=Zo  &  a^F^((T 

F  ( <7 ,  a  )  =  Zo  and  "1 
daF^(aya)&>0  j 


,Q)a  <  0  } 


(3.9) 

(3.10) 


For  the  above  two  dimensional  elasto-plasticity  problem ,  Bonnetier  [3]  has  proved  the  following 
theorem : 

Theorem  3.1  If  the  gauge  function  F{c,a)  is  admissible ,  i.e.  it  satisfies  all  the  three  assumptions 
AI  )'^  A3),  and  if  the  body  traction  functions  fi{x  ,y,t)  ,  :  =  1 , 2  ore  in  (0)  for  all  fixed  t  ,  and 

piecewise  analytic  in  t  for  (  z ,  y )  a.e.  on  0  .  Then  the  two  dimensional  elasto-plasticity  problem 
(3.2)  'v  (3.10)  has  a  unique  solution.  Moreover,  the  solution  is  also  piecewise  analytic  in  time  t  . 


8 


4  Semi-Discrete  Approximation  of  the  Quasi-Static  Problem 


Consider  first  a  special  case  where  Q  C 
is  a  bounded  domain  which  can  be  parti¬ 
tioned  into  rectangles  {Tn}  of  size  h  and  the 
aspect  ratios  satisfy  ao  <  a  <  ai  where  oq 
and  tti  are  independent  of  the  element  size 
h  .  Let  Tk  be  such  a  partition ,  A/*  be  the 
set  of  nodes  of  Q\d(l ,  and  {  }  i  »  €  ^ 

be  the  piecewise  bilinear  basis  functions  such 
that  ^  =  1  at  the  node  n  and  =  0 
at  other  nodes. 


Figure  4.1  A  domain  il  which  can  be 
partitioned  into  several  rectangles 


Denote  by  Qifi{Slh.)  the  set  of  contin¬ 
uous  piecewise  bilinear  functions  defined  over  Q  such  that  they  equal  zero  on  the  boundary  d  il ,  t.e. 


■{ 


h{x  ,y)\T„  =  ax  +  by  +  ex  y  +  d 
h(x,y)\Ba  =  0 


) 


Thus  any  function  H  (x  ,y)  in  Qijq  (D*)  can  be  written  in  the  form : 


n€V 

Meanwhile  denote  by  Qi  (fl^)  the  set  of  all  piecewise  bilinear  functions  defined  over  ,  i.e. 

Qi  (fill)  =  I  I  h(z,y)j^  =ai  +  6y-|-ciy-l-d|  . 

Notice  that  the  functions  in  Qi  (il/,)  arc  not  necessarily  continuous  over  the  domain  fi  .  Let  Af  be 
the  total  degrees  of  freedom  of  the  finite  element  space  Qi  {ilh) ,  and  {  }  ,  n  ^  Af  be  the  piecewise 

bilinear  basis  functions  of  the  space  (Dh)  .  So  any  function  H  {x  ,y)  in  (D^)  can  be  expressed 
in  the  form : 


n{x,y)=  51  Hn^n{x,y). 
n€V 

For  piecewise  bilinear  functions ,  we  can  use  four  Gaussian  points  on  each  rectangular  element .  Then  by 
using  the  interpolation  method  ,  any  function  in  the  space  Qi  (D/i)  will  be  uniquely  determined  by  its 
values  at  these  four  Gaussian  points . 

Now  we  can  formulate  our  semi-discrete  plasticity  problem  as  follows . 
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To  seek  the  solution  (Cf,<7,a)  ,  such  that  for  any  fixed  t  ,  (/=  (Ui,  l/j)^  €  (Qi,o(iih)f  > 
<7  as  (<7n.ffw,<yw)^  €  (^i(flA))  ,  a  a=  (o,  ,a3,...,ot„)^  €  (Qi(ft*))”*  ,  and  they  are  in  the 
form 

With  f/..„€C^‘*(o,r) 

II 

=  E  <^O.n(0^n(*) 

with  ffii.*  €C®'**(0,7’) 

II 

Oi  =  «i.n(0^n(*) 

n€A^ 

A  _... 

with  a..„€Ci*‘*(0,r) 

II 

•• 

• 

3 

where  (  0,T)  is  the  set  of  functions ,  continuous  on  [0 ,  T] ,  and  have  a  uniformly  bounded  right 
derivative  of  any  order  at  each  point  on  [  0 ,  T ) .  The  stress  functions  satisfy  the  equilibrium  equations: 


V  n  6 


(4.1) 


Here  the  body  traction  /  =  ( /i ,  /j )  is  piecewise  analytic  in  time  and  for  each  t  €  [0,r] ,  /(*,<)  € 
X’(n)  .  On  each  element ,  the  corresponding  constitutive  equations  are  given  by : 


d{xk)  =  Of  (it) 

a(**)  =  0 


if 


{  F{cr{xk),a{xk))  <  Zo 


or 


d„F^ixk)aixk)  <  0 


(4.2) 


Q  rX/ •\*/**\  *  ind 


(4.3) 


where  it  {k  =  1,2, 3, 4)  €  G  and  G  is  the  set  of  Gaussian  points  of  all  the  rectangles  in  the  dommn 
Q  ,  and 


The  initial  conditions  are : 


O'(it) 


Dd^Fd^F'^D 
daF^daF  +  d^F^Dd„F 


(4.4) 


[/(0,x)  =  0  ,  <7(0,  a:)  =  0  ,  o(0,a:)  =  0  .  (4.5) 

or  course  here  again  we  assume  that  the  traction  functions  satisfy  the  compatibility  conditions  : 

/i(0,a:)  =  /2(0,x)  =  0  . 


To  simplify  the  statement  of  our  proof,  in  the  following  context  we  will  assume  that  the  body 
traction  functions  are  analytic  in  time ,  and  the  gauge  function  F  ( <7 ,  a )  is  also  analytic  in  <7  and 
i»  .  ( Similar  results  still  hold  if  /  is  only  piecewise  analytic  in  time ,  and  F  ( er ,  o )  is  only  piecewise 
analytic  in  <7  and  a  .) 
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5  Existence  and  Uniqueness  of  the  Solution  of  the  Semi-discrete 
Problem 


Note  that  U  —  {U\  ,1/7)  is  chosen  to  be  piecewise  bilinear .  From  the  relation  between  strain 
and  displacement  (3.1),  we  know  that  the  strain  functions  as  well  as  the  strain  rate  functions  will  be 
piecewise  linear .  So  now  inst  of  the  whole  element ,  we  can  define  a  Gaussian  point  Xk  ^  G  of  a. 
rectangular  element  to  be  elastic  at  to  (resp.  plastic)  if  it  satisfies  (4.2)  (resp.  (4.3) )  at  the  point  xk  . 
Therefore  from  their  values  af  the  Gaussian  points  we  can  define  the  functions  tr  ,  a  ,  a  and  d  to  be 
some  piecewise  bilinear  functions . 

Consequently  we  can  now  define  the  set  of  elastic  Gaussian  points  and  the  set  of  plastic  Gaussian 
points  as  follows: 


x£G 


F(a,a) (x)  <  Zq  or 

F  (0,  q)  (x)  =  Zq  and  d„  F"^  D 


«(£/)(*)  <0} 


V  =  V{V)  =  {x€G  I  F(a,a)ix)=ZoajiAdaF'^D(iU)ix)>0}  ,  (5.2) 

where  e{U)  =  ~  ~  ^  ^  Gaussian  points  in  the  dommn  Cl . 

Using  the  constitutive  equations  given  by  (4.2)  and  (4.3) ,  we  can  easily  verify  the  following : 

Proposition  5.1  The  elastic  and  the  plastic  Gaussian  point  sets  ,  £  and  V  ,  defined  by  (5.1)  and 
(5.2)  can  also  be  expressed  in  terms  of  the  stress  tensor:  W  U  €  (Bo  (Cl)  )* , 


l(U)  =  I 


X  €  G 


F(<rj  a)(x)  <  Zo  or 

F(<r,a)(x)=:Zo  &  daF^&(U)(x) 


J 


V(U)=^xeG  I  F(o,a)(x)=-.Zo  and  d^F^&(U)(x)  >  oj 
where  the  stress  tensor  (r(U)(x)  is  defined  by  the  constitutive  equations  (4.2)  and  (4-3). 


Let  ( •  ,  • )  be  the  standard  vector  inner  product  defined  on  R®  .  For  any  given  displacement 
function  U  (x  ,i)  and  any  Gaussian  point  x  ,  by  using  (3.1)  we  can  easily  find  the  values  of  the 
corresponding  strain  rate  at  that  point,  e(x,t)  .  Then  by  using  the  constitutive  equations  (4.2)  and 
(4.3)  we  can  find  the  values  of  the  corresponding  stress  rate  at  that  point,  cr(x,t)  .  Hence  at  each 
Gaussian  point  x  ,  we  can  defined  a  functional  G  (U  (x) )  such  that : 

G(U(x))  =  (&(U(x)),i(U(x))). 

For  this  functional ,  we  can  give  the  following  Lemma  which  is  important  in  the  proof  of  existence  and 
uniqueness  of  the  solution  for  the  semi-discrete  problem . 
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Lemma  5.1  Let  {a{xk)tdi{xk))  €  x  2/  be  determined  by  equation  (4.2)  and  (4.S).  Then  the 
functional  G(U{xk))  =  (^iU{xk)),€{U{xk)))  ,  defined  on  the  eet  (Qifi{Tn))^  ,  w  strictly  convex. 
Here  is  any  rectangular  element  of  the  partition  %  . 

Proof:  Let  £2,  V  e  (Oi^(rn))’,  tf  €  (0, 1)  and  VP  =  tft/ +  (1  -  .  We  define: 

fb  =  {i*€r«nG|/’(<r,a)  <  Zo) 

Pi  =  {*i€r„nG  I  F(a,a)  =  Zo,d„F^Di{U)>ll,daf^Diiy)>Q} 

P,=  {ik€rnn<5  I  F{a,a)  =  Zo,dcF^Di{U)<ll,d„F^Di{y)<o} 
P3={ifc€r„nG  I  F{a,OL)  =  ZQ,d^F‘^Di{V)>(i,dcF^DHy)<ll) 

7>4  =  {iib  eTnOG  I  f’(ff,a)  =  2o,  5,J’^2)c(t2)<0,  d„F^Df(V)>0} 

1 ) .  On  €q\JV2  three  functions ,  U  ,V  and  VP  ,  are  corresponding  to  the  elastic  state .  So  we  get 

G(IP)  =  (o(VP),€(W))  =  (Dc(W),c(VP)) 

<  9iDe{U),€iU))  +  {l-0){DiiV),eiV)) 

=  0{diU),e(U))  +  il-9)iaiV),e(V)) 

=  0G(U)4-  {l-9)GiV)  . 

2  ).  On  Vi  ,  all  three  functions ,  £/ ,  V  and  VP  ,  are  corresponding  to  the  plastic  state .  Hence  we  have 

G(W)  =  (c(W),iiW))  =  ((D-D‘)iiW),i{W)) 

<  9aD-D')iiU),€{U))  +  {l-0)aD-D')i{V),i{V)) 

=  (?(d(£/),c(£/))  +  (l-2>)(d(P),c(P)) 

=  »G(l/)  + (l-e)G(P)  , 

where  D '  is  defined  as  (4.4)  which  is  independent  of  the  values  of  U  and  V  ,  and  only  dependents  on 
the  values  of  <r  and  a  . 

3  ).  On  Vz  ,  we  have  the  case  where  U  is  corresponding  to  the  elastic  state ,  and  V  corresponding 
to  the  plastic  state .  Decompose  the  vector  e  =  (  en ,  €32  >  ^13  into  a  sum  of  two  vectors ,  one  parallel 
to  d(,F  (<7  (x/t) ,  a  (x;t))  1  ^.nd  the  other  orthogonal  to  5„F(<7(xit) ,  a(xfc) )  with  respect  to  the  scalar 
product  (•,•),  i.e. 


€  =  C-^  +  C  * 

(5.5) 

and 

(De-^,€*)  =  0  . 

(5.6) 

Obviously  we  have 

r  €-^(iP) 

=  9c^iU)  +  ii-e)c^(v) 

< 

1 

1  i'm 

=  e€''(t/)  + (l-«)c"(P) 

(5.7) 
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Caiel.  StippoM  d,F^Di(W)  >  0,thea 

=  {€(tV)  -  ^».F^ DHW)a,F,HW)) 

+  {i\W)-Md^F'^D€iW)d,F,i*iW)), 

where  M  =  ( daF'^daF  +  doF^  D  B^F  )  ( <r  (**) ,  a  (xk)  ) .  Hence  by  the  definition  of  c  *  ,  we  have 
(<7(W),€(W))  =  {€^(W),e^iW))  +  KidaF^Di{W))\  (5. 


where 


^  ■”  [  d„FT  D  d^F  daF^  D  dcF  +  doF^  daF  ] 


is  a  positive  number  independent  of  the  values  of  U  and  V  .  FVom  the  fact  that 

(a^F^DiiW)  =  ed^F"^  Di{U)  +  il-e)d„F‘^D€iV)>  0 


d^F^DiiV)  <  0 


we  get 


0  <  a^F^DeiW)  <  ea^F'^DeiU) 


On  the  other  hand,  by  the  convexity  of  the  scalar  product  ( * ,  * )  and  the  equation  (5.7),  we  have 

<  9(€^iU),€^{U))  +  il~0)(€^{V),€^{V))  . 

Hence  the  equation  (5.8)  becomes 

(»(H'),4(H'))  <  »(r^(P).4‘(i;))  +  (l-*){c^(V),e'^(V)) 

+  »jr(a,F’'X)e(J7))’  , 

or 

(d(ir),f(H'))  <  0  \{i^{U).i^{U))  +  Kid^F'^Diiuyf] 

^  ■‘(5  9) 

+  (l-<i)(c(V),6(V))  . 

Using  (5.8)  with  W  replaced  by  U  (since  &{U)  corresponding  to  the  plastic  state ),  we  finally  have 
( d  (ly) ,  e  (ly) )  <  0  ( d  (U) ,  e  ({f ) )  +  (1  -  e)  ( d  (V) ,  e  (F) )  . 


Case2.  Suppose  dt,F^Di(W)  <  0  ,  then 

(d(lF),e(lF))  =  (€(1F),€(VK)) 
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or 


FVom  the  fact  that 


{d^F^DiiW)  a  ed^F^D€iU)  + {l-9)d„F^D€iV) 
d„F^Di{U)  >  0 


wc  get 


(1  -  6)  d„F^  DiiV)<  daF"^  D€(W)  . 

Hence  (5.10)  can  be  rewritten  as 

<  (l-»)  (t^(V),i^(V))  + 

Replacing  W  by  U  in  (5.8) ,  wc  get 

{i^(U),i^iU))s(HV),nv))  . 

Replacing  IV  by  V  in  (5.10)  (since  &{V)  is  corresponding  to  the  elastic  state ),  we  get 


( j(V).t(V))  =  (c^(V).<^(V))  -f 


So  again  we  get 


(5.10) 


i&(W),i(W))  <  (l-e)(a(V),€iV))  +  9{cr{U),€{U)) 


4).  On  V4  ,  we  can  prove  the  same  result  by  changing  the  roles  of  U  and  V  in  3). 

fVom  the  above  discussion  we  see  that  the  functional  G{U  {xk) )  is  strictly  convex ,  and  this  com¬ 
pletes  our  proof  of  the  lemma.  □ 

Now  let  /  =  ( /i ,  /a  )  €  {L^  (0)  )*  be  the  incremental  body  traction,  we  can  define  a  functional 
a.s  ro)low.s: 

For  any 

ir  =  [  |-  ]  6  (Hi(n))*  , 

let 

F{<r,a,U)  a  1  /  ia(U),HU))dx~  f  {f,U)dx  (5.11) 

^  Jtt  Jn 

where  ff(U)  €  (Oi(ftA))^  is  defined  as  in  (4.2)  and  (4.3).  Then  by  using  the  result  of  Lemma 5.1 ,  we 
can  easily  prove  the  following  theorem . 
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Theorem  6.1  The  funetioruU  ^(<7 ,  a ,  •)  defined  by  (5.11)  on  (  (rt)  )’  is  strictly  convex  and  coercive 
over  the  space 

Proof:  Becanse  ^{U)  is  piecewise  bilinear,  and  i(U)  is  piecewise  linear  over  0  ,  for  any 

rectangular  element  Tn  €  7^  ,  we  have  that  {a(U),€{U))  €  Q2 (Tn) ,  where 

h(»,  y)  =oi  +  oax  +  oav +  «4*y  +  osx’  +  a6X*y 

+07 X y^ +  o«y*  + <19 V(x,y)€r„ 

For  the  four  points  Gaussian  quadrature,  it  is  exact  for  any  function  in  Q2  >  i.e.  for  any  Tn  €  Tk 
have 

/  (d(ff),£(lf))dx=  ^a»*(d(t^(xfc)),£(t/(xfc)))  ,  (5.12) 

Jbsl 


■I 


<?2(rH)=U(x,y) 


where  it  (k  =  1,2, 3,4)  are  the  Gaussian  points  of  the  rectangular  Tn  ,  and  u/k  (k  =  1,2, 3, 4)  are  the 
weighting  coefficients. 

fVom  LemmaS.l,  we  know  that  the  functional  G((/(£k))  =  (^(f^(xjb))>c(f^(xib)))  is  strictly  convex 
over  (Qifl(Tn))^  .  We  also  know  that  the  weighting  coefficients,  Uk  (k  =  1,2, 3, 4)  ,  are  all  positive. 
Therefore  for  any  U  ,V  €  ( <?i.o  (flfc)  )^  >  6  €  (0,1)  and  W  =:  +  (1  -  ®)  V  ,  we  have 


f  (d(W).c(W))dx=  ^u;*(o(W(x*)),€(W(x*))) 

=  ^u;fe(<7(df/(x*)  +  (l-d)V(xfc)),c(^£/(x*)  +  (l-^)V(xk))) 

*ssl 

<  ^a;fcd(d(i; (**)),  f(f^(xife)))  +  ^u;*(l-<>)(a(V(x*)),c(V(x*))) 

IIrsI  fcsl 


=  <?/  i&(U),iiU))dx+(l-S)  f  i&iV),iiV))dx  , 

JTn  JT„ 

for  any  T„  £  Ti,  .  Obviously  /  f  U  dx  ]s  a,  linear  functional  of  U  .  So  .F ( a ,  a ,  • )  is  strictly  convex 

Jn 

over  {Qt,o(i^h)f  . 

To  prove  the  coerciveness,  take  a  sequence  {ifn)SILi  €  (Qi,o(nfc))^,  such  that  ||  Un  11^1(0)  — ^  o® 
as  R  — *  00  .  Then  from  Proposition  2.1 ,  we  know  that  there  exists  a  constant  7' 

iaiU{xk)),HUixk)))  >  't'iiiUixk))MU{xk)))  k  =  1,2, 3,4, 

where  7'  is  independent  of  the  value  of  a,a,i  and  U  (xk) . 
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Heaet, 

nt.a.O.)  .  jf 


E  E"»(*(P.(i»)).<(^(i*)))-  f  U.6.)dx 

T«€Tfc  k«l 


£  liLUn)d: 

r«€Tfc  Iwl 


=  I  jf  («(«,).  i(P.))<i*  -  jf  (/. 


by  unng  the  well  known  Korn’s  Inequality  ( see ,  e.g.  [4]  [13]  [8] ).  Hence  as  ||  Un  ||//i(n) 
T{aya,Vn)  — *  00  ,  and  the  coerdveness  is  proved. 


00  ,  we  have 
□ 


Now  let  us  consider  the  following  minimization  problem 

min  .F(a,a,y)=  min  ^  I (a{V)yi{U))dx-  [  fUdx  (5.13) 

(rettk  OeHk  ^  Jo  Jo 

where  Hk  =  (Qi,o(^a))^  •  Then  from  the  previous  theorem,  we  know  that  for  any  (<r,a)  € 

{Qii^k))^  X  I  is  a  strictly  convex  and  coerdve  functional  of  U  €  Hk  . 

So  the  minimization  problem  (5.13)  must  have  a  unique  solution  U  £  Hk  •  Therefore  we  can  have  the 

following  corollary: 

Corollary  S.l  The  functional  ^(0,0,')  defined  in  (5.11)  has  a  unique  minimum  over  the  space 

iQiMW- 


By  showing  that  the  functional  /* ( <7 ,  a ,  * )  is  fVdchet  differentiable  on  the  space  ( Qi,o  (Ha)  )  ^  1 
we  then  have  the  following  theorem : 


Theorem  5.2  The  minimum  U  of  (5.13)  is  the  solution  of  the  following  problem: 

dit>n 


n  G  Af 


F{(T{xk),a(xk))  <  Zo 


or 


daF^{xk)oixk)  <  0 


with  the  constitutive  equations: 

a{xk)  =  Deixk) 

a{xk)  =  0 

fT{ik)  =  (D-  D){xk)iiik)  (  F{a{xk)yQixk))  =  Zo 

here  Xk  €  G  is  any  Gaussian  point  in  the  domain  U  ,  and  D'  {xk)  is  defined  as  in  (f.4). 


"  1 


(5.14) 


(5.15) 


(5.16) 
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PROor: 


(7  ,  ,  1/  >  0  and  W  =  U  +  .  Denote 

(^S»hi("))  -  (.(g).iw)] 


E“*( 

*•1  *- 


(diW,(xk)),HWAik)))  -  ciUiik)),HUiik))) 


2v 


-(i(»(f4)), <(*(*,)))]  . 


and 


We  will  discuss  the  fdlowing  live  different  cases: 

ik  €  €o  ={*€5  I  P(cr(*),0((»))  <  Zo} 


ik  €  pr 


■1 

i*  €  n  =  1 


x€G  I  f’(<r(*),ot(*))  =  Zb  , 
X  €  5  I  P(o(x),ot(x))  ss  Zb  , 


ik  € 


S  I  X  €  ^ 


I  F(o(x),a(x))=  Zo, 


■{ 


i*  €  7?;  «  ^  X  €  C?  I  /’(«r(x),o(x))  *  Zo  , 


>.j 

>0  J 

?«(£f)(x)<0  I 

'«(W'.,)(*)<0  J 

'€(t^)(x)>0  I 

«(W'..)(*)<0  J 

'<(f^)(»)<o  I 

c(W„)(x)>0  j 


ft,f^D£(y)(x)  >  0 

doF^  Di{U){x)<{i 
deF'^D 

d*F^D€(t^)(x)>0 
d,F^D 


d„F^D€{U)ix)<0 
d„F^D 


1).  On  the  set  £o  U  ^3  > 


or 


A.(x)  =  -|(Df(^),e(^))  . 


2).  Similarly  on  the  set  Vi  ,  we  have 
K{x) 


=  -!((- 


Dd„Fd„F’^D 
daF^d„F-\-d„F^Dd„F 


3).  On  the  set  ,  we  can  use  the  decomposition  (2.15),  (2.16), (2.17)  and  get 

Mi)  =  ^  iid„F'^Di{V)d,F ,€{U))) 


(5.17) 


(5.18) 
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-WV)-iid.f^DtiV)d.F,H*)) 

(/lAf ,  <(U))+  tiiia.F,  lit)) 

here  7  ss  d^F^Di^U) .  So  we  get 

F^m  the  definition ,  we  know  that  on  Pjf  we  have 

f  a^I^Di{V)  >  0 
\  d^F^D€(U)-k-ud^f^Diiifi)  <  0 
Hence  I7I  =  \d„F^Di{U)\  <  1/ 1  £ (^) | ,  and  (5.19)  yields 

h.(x)  <  f  <  r (^) , £(^)  >+ 1|  £ (^)  |»  ( > 

3).  On  the  set  V4  ,  we  have 

.  {HW.)-,id^F^DHWMFMW,))-{€{U)MU)) 

- ^ - - - 

Hm.im 

=  ■^d,F'^Di{W,){d„F,k(W,)) 

=  -  ^d,F'^Di{W,)(,d,F,i{JJ)-^vi{il>))  . 

From  the  dcRnition  of  set  ,  we  have 

J  d„F'^D({U)^-vd„F'^Di{il>)  >  0 
\  d„F'^Di{U)  <  0 

which  yields  |  d^F^ D  i  (U)  |  <  t'  |  daF'^D  £  (^)  | .  So  finally  we  get 

IMi)l  <  +  ^{d,F,i{U)  +  ui{<l,))^ 

if  4»U 

=  +  ~  ((a.F,£(l/)>V 

2v{^„FMV)){^„F,m)■^u^^„FM<t>)?) 


(5.19) 


(5.20) 


(5.21) 


18 


Therefore  using  (547), (5.18), (5.20), (5.21)  and  the  uniform  boundedness  of  the  gradient  d^F ,  we  can 
get  that 

\hy{ik)\  <  i^C(£(^(fs)),€(^(x*)))  €  if*  and  ifc=  1,2,3,4 

for  s(une  C  >  0 . 


Hence  finally  we  get 


Fia,a,W^)-F(a,a,U) 


Wfc  (**) 


kml 


<  2  vCw*(<(^(*k)),£(^(**))> 


kml 


=  I/Cjf^(l?€(^(x)),£(^(*)))dx, 


which  means  F(a,a,-)  is  differentiable  on  (Qi^(flk))^.  Differentiating  F  with  respect  to  C/,  we  get 
that  the  minimum  U  of  (5.13)  satisfies  the  system  (5.14),  (5.15)  and  (5.16).  This  completes  the  proof  of 
the  theorem .  □ 


By  using  the  results  of  Corollary  5.1  and  Theorem  5.2,  we  can  now  prove  the  following  theorem 
about  existence  and  uniqueness  of  the  solution  of  the  semi-discrete  problem. 

Theorem  5.3  If  the  body  traction  function  f  (x,f)  is  piecewise  analytic  in  t  ,  and  for  every  fixed 
t  €  [0,T] ,  f(x,t)  €  X^(D) ,  then  there  is  a  unique  sokiion  for  the  ODE  system  (4-1)  ,  (4-i)  ,  (4’S) 
and  (4-4)  satisfying  the  initial  condition  (4-5),  and  moreover  the  solution  is  piecewise  analytic  in  time 
t  . 


Proof  :  We  will  prove  this  theorem  in  4  steps .  For  the  sake  simplicity ,  we  will  assume  that  f  ( x ,  f ) 
is  analytic  in  t  and  F(o,q)  is  analytic  in  o  and  a  -  But  the  same  result  still  holds  when  f  and 
F  are  piecewise  analytic . 

Stepl .  Here  we  will  assume  that  at  the  beginning  the  material  is  in  its  virgin  state.  So  at  all  the 
Gaussian  points  the  material  is  elastic  until  a  time  t  =  to  when  F(o,a) (it ,to)  =  F(o,  0) (it , to)  = 

A  A 

Zo  at  some  Gaussian  points  of  some  elements.  Denote  by  Gj  the  set  of  those  Gaussian  points  and  G 
the  set  of  all  the  Gaussian  points  of  Qt  .  Then  for  t  €  [0,to)  ,  all  the  points  are  elastic,  and  we  just 
have  a  cla.<;sical  linear  cla.sticity  problem.  So  the  problem  wiU  have  a  unique  analytic  solution  over  the 
time  interval  [0,to). 

A  ^ 

Step 2.  As  we  will  see  in  the  following  proof,  the  points  of  the  set  G\Gi  will  remain  elastic  after 
to  .  But  those  of  Gi  may  or  may  not  yield  after  to  .  To  determine  this  we  have  to  check  the  sign  of 
d„F^{a ,a)a(xk  ,to  +  0) .  Consider  an  auxiliary  equation  system  such  that  (a^  ,U^)  €  (Qi(D/i))^  x 
( QiA^k)  )^  would  solve  at  to  +  0  . 
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X<"'”  “  X*^'****"’’**’ 


n  €  Af 


(5.22) 


with  the  amititvtive  eqve^mt: 


<y*(i*)  s=  D€^(ik) 


ff^iXk)  =  /?€*(**) 


(ifc)  = 


Vi*  6  G\di 


Vi*  € 
Vi*  € 


(5.23) 


(5.24) 


where 


i>l  =  {x€di  I  d*F^(«o)X>€(l^*)(i)  <  O} 

^;  =  {i€<?i  I  d^F^ito)D€iU')(x)  >  O] 

and 

n‘  -  dV,  m  -  Dd.F».f’D  .  . , 

FVom  C<»oIlary5.1  and  Theorem 5.2,  we  know  that  there  exists  a  unique  solution  of  the 

equation  system  (5.22)  ~  (5.24). 

Step  3 .  Now  we  can  define  the  sets 

^={i€Gi  I  d,,F^(i,to)<^M»)  >  O} 

^  =  {i€G,  I  a<rF^(i,<o)‘^*(i)  <  O} 
and  consider  the  following  initial  value  problem  starting  at  t  =  to’. 


^(a„,^)  +  (<^«.^)  =  l{f2{t)An) 
f  d(i*)  =  X?£(i*) 

{  Vi*eG\Gj 

I  d(i*)  =  0 


d(i*)  =  {D-D')iixk) 

A(t  \  i^k) a jxk)  ^  . 


Vi*  €  ^ 


(5.25) 


(5.26) 


(5.27) 
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with 

i(f,)  =  <(tf(iO) 

,  Da.F^ik)d.f^(it)D 

'**'  ■  a.f’-  (i»)9./'(i») + i>«,f  (i») 

and  the  initial  conditions : 


(y,o,o)(»*.<o)  =  «)(**. <o-0)  .  (5.28) 

Now  we  can  show  that  the  initial  value  problem  (5.25)''«(5.28)  has  a  unique  analytic  solution  {U,a^a) 
on  the  interval  [to  +  ^)  with  some  j  >  0 ,  and  moreover 

(U,^)k+o  =  (U^oM 

where  ( )  is  the  solution  of  the  problem  (5.22)  '<-(5.24). 

Indeed ,  substituting  (5.26)  and  (5.27)  into  (5.25) ,  because  of  the  assumption  that  the  gauge  function 
F{tr  f  a)  is  an  analytic  function  of  c  and  a  ,  we  will  get  an  equation  in  U  with  analytic  coefficients. 
On  the  other  hand  from  Proposition 2.1 ,  we  know  that  D  and  D  —  D'  are  positive  definite  matrices. 
Therefore  there  exists  a  unique  solution  V  which  is  analytic  on  a  interval  [  to » fo  +  j)  for  some  j  >  0  . 
For  the  stress  and  the  hardening  parameter  functions  ( <r ,  a  ) ,  again  because  of  the  analytic  assumption 
of  the  gauge  function  F(<r  ,o)  ,  we  can  write  (5.26)  and  (5.27)  as  a  system  of  ODE’s  in  the  form 

~  s=  ^(<T,a,t)  (5.29) 

[  o 

where  A{a,a,t)  is  an  analytic  function  of  <r ,  a  and  t .  Therefore  (5.29)  would  have  a  unique  analytic 
scdution  over  the  interval  [  to  •  to  +  ^  )  satisfying  the  initial  conditions  (5.28) . 


The  solution  of  problem  (5.22) 'v  (5.24),  ( t/*  ) ,  satisfies  the  following  problem  at  to  +  0 


n  €  -V  (5.30) 


where 


<T*(x*)  =  D(^{Xk) 


Vie  G\G\ 


(5.31) 


[  a^ixk)  =  (£>-/?')€»  (ifc)  Vi  €  G? 

with  17'  e  (Qi,o(flfc)  )^  •  By  virtue  of  Theorem 5.2,  this  problem  has  a  unique  solution.  On  the  other 
hand,  {U ,  a)  ,  the  derivative  of  the  solution  of  problem  (5.25)  ~  (5.28)  with  respect  to  time  t  also 
satisfies  (5.30)  ~  (5.31)  at  to  +  0  ,  therefore  we  have  {U ,  a)  |<o+o  =  ( t/' ,  <t'  ) . 


Step  4 .  In  this  step  we  will  show  that  there  exists  a  constant  >  0  ,  such  that  at  any  Gaussian  point 
Xk  ,  if  it  is  elastic  (or  plastic  )  at  the  time  t  =  to  +  0  ,  then  it  will  still  be  elastic  ( or  plastic )  over  the 
time  interval  [  to ,  to  +  ^ ) . 
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a) .  Fbr  any  Gauisian  point  it  €  G\Gi  =  |  £  €  G  |  F(a^  ,a^)(.T)  <  ZoJ,  becaune  the  solution 

(17,9,0)  is  analytic  over  [ to •  <o +  £),  to  we  can  always  find  a  constant  0  <  <  d  such  that 

f’(<r,a)(xA,t)  <  %  V£*€G\Gi  over  the  interval  [to>to  +  ^i) 

which  iinpUes  that  the  Gaussian  pdnts  in  ^Gi  will  remain  in  elastic  state  over  the  time  interval 
{<o.<o  +  ^) . 

b) .  For  any  Gaussian  point  i*  €  G  such  that  F(a^  ,Q*)(ifc)  =  2b  and  <  0  , 

we  have 

•  .1  /4+*'  _  _ 

F*((r,a)(i*,to  +  i')-/’’(<T,a)(*fc,to)=2^  Fia ,a){d^F^&  ■¥  d^F^a)di  . 

So  if  If  <  6  is  small  enough,  we  have  d^F^ a{xkti)  <  0  and  d(x;k*0  =  0  for  any  t  €  [tot<o  +  t') . 
(from  the  analytic  assumption  of  the  gauge  function  F(<r,a)  and  the  analyticity  of  the  sedation 
functions  (f/,<r,a)  and  the  fact  that  9(xfc,to  +  0)  s  <r*(xik)  )  Therefore 

F^{(r,a){xk,to  +  v)- F^{<r,a)(ik,io)  =  \  f  F{a ,a)dgF'^ adt  <  0  , 

^  Jio 
i.e. 

F^{a,a){xk,to  +  tf)  <  F^ (c,a)ixk,to)  =  . 

Hence  for  t  €  [to«<o  +  we  have  F(9,a)(xjk,t)  <  2o  ,  which  means  that  the  elastic  Gaussian 

A  m 

points  of  Gj  remain  elastic  over  the  time  interval  [to  ,to  + 1^2)  • 

c) .  For  any  xt  €  G|  ,  we  also  have : 

A  n  1  r*o+«'  _  _ 

F^i<r,Q)ixk,io  +  *')- F^{a,a){xk,to)  =  -  f  F{a ,a){d„F^a  +  daF^a)dt  . 

Since  on  Gj  we  have  d„F^ o{xk,io  +  0)  >  0  ,  once  more  by  the  analyticity  of  the  gauge  function 
F(ff,a)  and  the  continuity  of  the  solution  9,0,9  and  d  ,  we  know  that  on  [^Ot^o  +  Wt^or  if  small 
enough,  we  have  d„F^ a  >  0  ,  and  moreover  d  cannot  vanish ,  so  F{a ,a){d„F^ 0  daF"^ a)  =  0. 
Hence  we  have  F(9,a)(x*,<)  =  F (9 , a) (x* , to)  =  2o  ,  V  t  €  [to,to  +  ^3)  >  which  means  that  the 
plastic  Gaussian  points  of  G\  remain  plastic. 

So  if  we  denote  by  Gj  the  set: 

G2  =  {x€G,  I  5^F^(x)9Mi)  =  0}  , 

and  if  Gj  is  empty ,  then  from  the  above  discussion  wc  can  see  that  the  elastic  Gaussian  point  set  and 
the  plastic  Gaussian  point  set  determined  by  the  solution  (  ,9^  ,a^ )  of  problem  (5.22)  ~  (5.24)  will 

remain  unchanged  for  a  certain  period  of  time  [  to  ?  ^0  +  ^  )  .  Thus  whenever  the  elastic  and  the  plastic 
Gaussian  point  sets  were  determined  at  the  time  level  t  =  to  +  0  ,  we  can  always  solve  the  initial  value 
problem  (5.14)~(5.16)  starting  at  t  =  to  and  forward  to  the  time  level  to  +  ^  without  violating  the 
admissible  conditions  F  (9 ,  a)  (x;^  <  t)  <  2o  • 
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Remark  5.1  If  is  not  empty,  we  can  show  that  the  choice  for  the  Gaussian  points  of  G2  will 
not  change  the  state  of  the  points  of  G\G2  •  By  considering  the  sign  of  higher  order  derivatives  of  the 
function  ,  to)  De  (x/b  *  to  +  0)  ,  we  can  also  show  that  there  exists  a  positive  number  v  ,  such 

that  for  t  €  [totto  +  t  the  elastic  Gaussian  points  remain  elastic  and  the  plastic  Gaussian  points 
remain  plastic.  The  details  of  the  proof  can  be  obtmned  by  slightly  modifying  the  proof  in  Bonnetier’s 
thesis  [3] . 

In  conclusion,  we  can  always  determine  the  state  of  each  Gaussian  point  of  the  set  G  for  a  certain 
period  of  time  after  to  by  only  knowing  the  history  of  the  solution  before  to  and  the  information  about 
the  traction  function  /(x  ,1©  +  0) .  Therefore  the  problem  (5.25)'^  (5.28)  always  has  a  unique  analytic 
solution  over  some  time  interval  [to<^+^)  with  the  sets  Gj  and  G\Gi  being  modified  a  little  bit  if  the 
set  Gj  is  not  empty.  Thus  we  can  always  solve  an  initial  value  problem  starting  at  t  =  t,-  and  forward 
to  another  time  level  t  =  tj+j  when  we  have  some  Gaussian  points  switch  their  elastic  or  plastic  state . 
Then  we  wiU  use  the  same  method  to  determine  the  elastic  or  plastic  state  for  each  Gaussian  point  for 
the  time  t  =  t,>i  +  0  ,  and  solve  another  initial  value  problem  starting  at  t  =  ii+i  .  The  similar  results 
still  hold  for  the  subsequent  yieldings.  Hence  finaUy  we  will  get  a  unique  solution ,  piecewise  analytic 
in  time ,  of  the  system  (4.1) ~  (4.4)  over  the  time  interval  [  0 ,T  ]  satisfying  the  initial  value  condition 
(4.5).  Thi.s  completes  Onr  proof  of  the  theorem.  D 

6  Energy  Estimates 

To  prove  the  convergence  of  the  solution  for  the  semi-discrete  problem ,  we  need  the  uniform  bound¬ 
edness  of  the  solution  .  By  using  the  admissible  assumptions ,  we  can  have  the  following  energy  estimate. 

Theorem  6.1  There  is  a  positive  constant  C  independent  of  the  finite  element  mesh,  such  that 

ll^llw^(n) »  l|f|lL»(n) .  Il^llt2(n)»  ll«lli,»(n)  ^  C’||/||i2(n) 

o.c.  t  €  (0,r) . 

Proof:  Let  {lJ,o,a)  be  a  solution  of  (4.1) (4.4)  with  the  initial  condition  (4.5)  in  the 

interval  [0,7)  .  Then  at  any  Gaussian  point  x*  ,  from  Proposition  2.1 ,  we  know  that  there  exists  a 
constant  c  independent  of  the  element  mesh  such  that 

c  (c(xfc),e(xfc))  <  (^ (**).«(**))  'i  Xk  e  G  . 

So 

cl  ie,e)dx  =  c'^Uk{Hxk),i{xk))<Yiuk{^(xk),i{xk))=  /  {&,e)dx  . 

*=1 

By  virtue  of  (4.1) ,  we  get 

c  f  ii,€)dx  <  f  {a,()dx  =  f  {f,U)dx  .  (6.1) 

JQ  Jo  Jo 
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Unng  Korn’s  inequality,  we  have 


II  ^  ^  llL»(n)  ^  ^  II  /llt>(n)  II  ^  lli,>(n)  • 

So  we  finally  get 

II  ^  ll«J(n)  ^  II  /llL»(n)  •  (6-2) 

FVom  (6.1)  and  (6.2),  we  have 

IUIlL>(n)  ^  II  /  llL»(n)  II  ^  lli*(n)  <  ^ill/lli»(o)  • 

For  the  boundedness  of  the  stress ,  we  have 

i  (**)  =  Z)“*  a  (i*)  if  is  elastic  , 

or 

«(**)  =  —  D')~^  a{xk)  if  Xk  is  plastic  , 

where  both  D~^  and  (D  —  D‘)~^  are  positive  definite ,  and  their  smallest  eigenvalue  is  bounded  below 
by  a  positive  constant  C]  independent  of  the  value  of  t,<r,a  and  the  element  mesh.  So  we  can  use  the 
same  argument  as  the  one  we  used  before ,  and  get 

C2  /  (o,o)d*  <  /  (f,U)dx  . 

Ja  Ja 

So 

II  ^  111,2(0)  ~  ^  II  llz,2(0)  II  U  lltaro)  <  C2  II  / 11^8(0)  ,  (6.3) 

by  using  the  inequality  (6.2) . 

For  the  boundedness  of  the  hardening  parameter ,  ||  a  ||£^(n) ,  we  can  use  the  fact  that 
d  (xk)  =  0  if  z;t  is  elastic  , 


d^r^(xk)a(xk)  « 


f.  /-  \]2  ^  O  (^/r)  a  r,/--  \ 

[(*)]< 


Vxk  €G  . 


Using  the  admissible  condition  7  <  |  | ,  |  dc,F  |  <  F  ,  we  then  have 


lo(*fc)]’  <  C  [d(iifc)f  , 


that  yields 


f  ^  f 

/  (d,Q)di  =  ^u;ifc(d(ifc),d(ifc))  <  C^u;jfc((T(zfc),d(iit))  =  C  {&,a)dx  . 

fc=i  fc=i 


By  virtue  of  (6.3) ,  it  follows  that 


/  (d,d)di  <C  f  (d,<7)di  <C  [  (/,/)d: 
Jn  Jn  Jn 


Hence  ||  d  ||£a(n)  <  C  ||  /  ||i,a(n)  and  the  proof  is  completed. 
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7  An  Equivalent  Weak  Formulation 


The  following  weak  formulation  is  the  key  point  of  our  proof  for  the  convergence  of  the  semi-discrete 
problem  solutions .  Based  on  this  formulation,  we  will  be  able  to  identify  the  equations  satisfied  by  the 
limit  of  the  semi-discrete  problem  solutions. 

Denote  the  set 

r  r*  =  2  rn(t)  .  r,(t)  €  01(1)  F(n,»'h)  (ik)  <  Zo ' 

^  I  /  \  neJf  with 

V  ngV 

Consider  then  the  functions  of  the  following  forms 


t/fc  =  E  ‘'“(OA.W  i=1.2 

n€V 

(7.1) 

ll€Ar 

<^h,j  =  ^  <t5(0^„(x)  *,j  =  1,2 

n€.V 

Ohi  =  £  «?(0^n(a:)  t  =  l,2,...,rii 

where  f/,"  (t) ,  o,"  (<) ,  oj*  (t)  €  C°’*®  (/) .  we  have  the  following  theorem. 

Theorem  7.1  The  semi  -  discrete  initial  value  problem  (4-1)'^  (4'^)  **  equivalent  to  the  following  prob- 


Seek  {Uh,ci,,ah)  of  the  form  (1.1)  such  that  for  all  t  £  I , 

(  f  ,  .  d<t>n  ,  .  X  .  [  t  i  J. 


L  ^ ^  =  /n  ( 


n 


where 


(<fh,Oh)  €  Kh 


Uke(Hl(il))^  Vt€/ 


J^{h-  Cah,Th-  ah)dx+J^(-ak,vh  -  oa  )di  <  0  V  (  ,  t/* )  € 


i«th  C  =  D~^  ,  and  the  initial  conditions 

[  £//.(0,*)  =  0 
^  <7fc(0,x)  =  0 
I  Ofc(0,i) =  0 


(7.5) 


Proof:  i)FiT8t.to8how  that  the  soliition  {V  ^a^a)  ofthe  system  (4.1)'^  (4.5)  satisfies  (7.2)  ~  (7.5). 

Indeed,  we  have  for  all  t  €  /  , 

F(<y(i*),o(z*))  <  Zo  'iik  ^  G  , 


which  means  ( tr ,  a  )  €  Kk-  Now  suppose  at  the  point  Xk  the  material  is  elastic,  then  from  (4.2) ,  we 
have 

f  i{xk)  «  C&ixk) 

\  d(xfc)  =  0 

Obviously ,  in  this  case ,  we  have 


(€(*fc)  -  Cd(i*).^(*fc)  -  <^(**))  +  (-dt(*fc),i^(i*)  -  a(z*))  <  0  .  (7.6) 


If  at  the  point  Xk  the  material  is  plastic,  then  from  (4.3) ,  we  know  that  there  exists  a  constant  A  >  0  , 
such  that 


f  «”(**)  =  Ad»f’(ff(i/k),a(xfc)) 


(7.7) 


(  -d(x*)  =  XdaF(tr(xk),a(xk)) 

where  (x*)  =  € (xk)  —  Co (xk)  .  On  the  other  hand,  since  F{c,a)  is  convex  and  analytic ,  we  have 
(  See,  e.g.  [18]  ) 


F{t^v)  -  /'(<7,a)  >  doF{ff {t  -  <t)  +  fieiP (a , a )^ ( i/  -  a)  . 


So  if  at  the  point  Xk  the  material  is  plastic  and  (r,i/)  g  Kk  ,'we  must  have 

F{T{xk),v{xk))-  F{a{xk),a(xk))-  F(T{xk),t'ixk))-  Zo  <  Zq-  Zo  =  0  . 


Hence 


d„F{<r  (xfc)  ,a(xfc)  )’^(r  (x*)  -  a  (x*) )  +  doF(a(xk) ,  ar(i*)  )^(i'  (it)  -  o  (xjt) )  <  0  . 


From  (7.7),  we  can  see  that  (7.6)  still  holds.  Therefore  on  any  T„  ^  Th  ,  we  have 

/  (i  —  C  a  ,T  -  er)dx  +  f  -  a)dx 

JTn  JTn 

k=i 


So 


f  ((-C&,T-er)(lx  +  f  (  -a  -  a)dx  <  0  V(r,i/)€  Kh 

Jn  JQ 


ii )  Second  to  show  that  the  solution  of  the  system  (7.2)  ^  (7.5)  is  unique.  Now  suppose  we  have  two 
solutions,  say  (tfi,ffi,ai)  and  (f/j  ,<^2  *03) « then  by  the  definition  we  have 


{<ri,ai)eKH  Vt€/ 
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( <73,03)  €  Kh 


By  virtue  of  (7.4) ,  it  fdlows  that 

/  ( ( ”  <2 )  —  C  ( Oi  —  03 ) ,  <73  —  oi )  <fa  +  ( ~( <ii  —  03 ) » oa  —  oj )  rf*  <  0  , 
which  implies 

^(«j-<3.<^3-<yi)rfx  +  [||V^(03-<^i)||’]+|^[||a3-a,||’]  <0  . 

FVom  (7.2),  we  can  show  that 

/  (€1  -C3,<73-<7i)<l®  =  0  . 

JO 


Therefore ,  we  have 


^  [||>/7(»j-<t,)||’  +  l|o,-o.||’]  <  0 


We  already  know  that  ( 17i ,  oj ,  oj )  and  ( f/3 « <^3 « 03  )  satisfy  the  same  initial  conditions  (7.5) ,  so  we 
get 


which  means 


*]  A  <  0  , 

i  "I  =<’> 

1  fl.C. 

1  «!  =  03 

fix/  . 

(7.8) 

/  <7,  =5  03 

1  A  A  ‘**«* 

1  O]  s  O3 

Slxl  . 

(7.9) 

It  follows  that 


For  the  strain  and  <iisplacement  function ,  we  know  that  at  any  Gaussian  point  Xk  ,  either 
€  (*fe)  =  D~^  a  (ifc)  if  ifc  is  elastic 


or 


Dd„  F  (<7 ,  a)d„F  (<7 ,  a)’^D 


doF{a,  ay^d^F  (<7 ,  o)  +  (<7 ,  ayDd„  F  (<7 

if  ik  is  plastic.  So  by  (7.8)  and  (7.9),  we  have 


.0)) 


-1 


^(®fc) 


•»=** 


<1  =  «3  . 


and  hence 


Ux=U2  . 

iii )  Finally  to  show  that  the  solution  of  the  system  (7.2) (7.5)  also  satisfies  (4.1)  ~  (4.5).  Indeed,  we 
have  for  all  <  €  /  , 


a)</z  +  I  {~a,u  -  a)dx  < 
JQ 


V(r,»/)€7u 
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So  for  any  GoaaiiMi  point  ik  t  take  r(f)  s  a{i)  ,  u{i)  =  o(£)  for  all  i  £  G  and  x  ^  it ,  and 
F(r(ia),i^(ijt))  <  Zb  for  £  s  .  Then  obviously  (r  ,t/)  €  .  and  so 


(i(£fc)  -  Cd(it), T(ik)  -  o(xfc))  +  (-d(ifc),i/(*fc)  -  o(ifc))  <  0  . 


(7.10) 


If  (£k»0  aatisfies  I'(ff(£ktt),oi(£kt*))  <  Zq  ,  by  the  continuity  of  the  gauge  function  F(a,a), 
we  can  conclude  from  (7.10)  that 


j  i(ik,t)^C&(xk,t 

[  d(£fc,t)=  0 


a{xk,t)  =  JDf(£*,t) 

< 

d(£fc,t)  s  0 


(7.11) 


which  is  just  (4.2). 


If  (£k ,  t)  satisfies  f  ( o  (£k « 1)  >  a  (xk ,  t)  )  —  Zq  ,  then  from  the  uniform  boundedness ,  |  a  (xk)  |  < 
C||/||£,3(0)  and  |d(£k)|  <  C  ||  /  ||£,3(o)  ,  we  know  that  a  and  a  are  absolutely  continuous  with 
respect  to  time  i.  So  if  we  define  g(t)  s  f'(<r(£k,t),a(xk,t) ),  because  F{a,a)  is  analytic  with 
respect  to  a  and  a ,  p(t)  is  differentiable  a.e.  on  [0,T] .  Hence 


lim 

g(t  +  h)  -  g(t) 

<  0 

because 

A — >0 

A 

lim 
h  *0 

g(i)  -  g(t-h) 
A 

>  0 

because 

and  so  we  get  (t)  —  0  ,  which  means 
at 

d»F(o(£k,l),a(£k,t))^d(£k,t)  +  5«F(o(£k,t),oi(£k,0)^d(*k,0  *  0 

Inequality  (7.10)  implies  that  there  exists  A  (xk ,f )  >  0  satisfying 


(7.12) 


£(xfc,0-Cd(xk,t) 

-d(xfc,t) 

From  (7.12)  we  get 


-d(xk,0 


dffF(a(£fc,l),a(ifc,t)) 

doF{<T{xk,t),a{xk,t)) 


A(xfc,0  = 

So  if  A  (xk  ,0  >  0  ,  we  have 


0<,F(ff(xk,t),a(xk,t))^g(xk,f) 


doF  (<r(xk,i),  o(xk ,  t )  )^daF  (a(xfc ,  t ) ,  o(xfc ,  t ) ) 


Hence 


d<rF((T(ik,i),o(xt,t))^&(xt,i)  >  0  . 


< 

.  /.  xx  dcF(<r,a)'^&  o 
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Solvlag  for  b  {ii, ,  t)  in  the  first  eqnntion,  we  get 


i,(A  —  ( n  _ Dd^FdsF^ D _ .  .  4\'tA  *\ 

&{xk,t)  -  {D  ^  d^FTDd^F  ^ ^ ‘ ^ 

a(xL.  t\  —  *^)^*^(***^)  ft  F(A  i\ 

o(**.0  -  {dc.F^dc,F)(xk,i) 
which  is  just  nothing  else  but  (4.3). 

If  A(£;k,t)  s  0  ,  then  we  have 

In  this  case,  we  again  get  (4.2)  and  the  proof  is  completed . 


8  Convergence  Analysis  of  the  Semi-discrete  Problem 


Let  (Uk^ffinOih)  he  the  unique  solution  of  the  semi-discrete  problem  (7.2)  (7.5)  associated  with 
the  partition  7a  of  the  domain  Q  .  Then  we  have  (Uht<rk,oih)  satisfying; 


n  eM' 


/  ((€fc-C^^fc,r-<^fc)  +  (-afc,i'“aa)]d®dt  <  0  'i{T,v)€Kk 

JQxI 


Uh  €  (gi.o(«/.))"  c  (^J(n))* 


for  all  t  €  / 


and  the  initial  conditions 

Uh{0,x)  =  0 
-  <7a(0,i)  =  0 
Qa(0,x) =  0 

From  Theorem  6.1 ,  we  know  that 


^Uh,Vk,€k,Qh 
^  Uh,Vh,ih,0‘h 


are  uniformly  bounded 


VUH,Uk,h,aH  J  inl~(/,l2(ft))  . 

Therefore  we  can  extract  a  subsequence  (  See ,  e.g.  (15)  &  [16]  ) ,  such  that 


^  Uh,  Uh,  (Th,  Ok  ,  ih 


W,U,tj,a,t 


^Uh,Uh,ah,Oh,ih 


V  U ,  U ,  & ,  a ,  € 


weakly  *  in  L°°  (/ ,  L^(fl)  )  .  For  the  limit  functions  <t  ,  o ,  c  and  U  ,  we  have  the  following  theorem . 
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Theorem  8.1  The  limit  functions  ( 1/ ,  a ,  o  )  is  the  unique  solution  of  the  follotaing  problem 


v«€  cs°(n) 


o.e.  <  €  /  with  (cr.o)  €  iT  ,  where 


K^{(<r,a)  €  L^(I,l\il))\F(c,a)ix,t)<Zo  a.e.  ftx/}  , 


and  the  admissible  equation 


J  ^  (c  -  Co,r  -  <r )+ ( — a)j  dxdf  <  0  V(r,i/)€A’ 


U,V  e  ,  a,o.a,o  €  , 


and  the  initial  conditions 


t^(*,0)  =  0 
<r(*,0)  =  0 
o(*,0)  :=  0 


a.e.  on  il 


Proof  ;  We  will  prove  this  theorem  in  five  steps. 

Step  1.  To  show  that  the  hmit  functions  ( If  ,(r,a)  also  satisfy  the  initial  condition  (8.8).  When 
choosing  the  subsequence  from  iUh,Oh-,Oh)  we  can  also  assume  that 

Uk,Uk  — - — ‘  U ,U  weakly*  in  I®® (/ , JTq (fi) )  »  (®-®) 

and  Uh  - -  U  weakly  in  (fix/), 

which  implies  that  UhtV  €  AC  (/ ,  )  » (  See,  e.g.  (15]  Chapter  11  ) .  So  for  any  <f>  (x)  €  (fi)  * 

t  —  T 

we  can  take  if{x,t)=  “jT"  ^(*)  5  8®^ 

V^(x,0€i*(/,X’(ft))  , 

Since  U  {x  ,t)  €  AC  {I ,  L\Q) )  ,  t/  (x ,  0)  is  well  defined  and 

U(x,0)<i>ix)  =  -  U  {x  dt  -  U(x,t)tlf(x,t)dt 

a.e.  X  in  f2  .  So  we  get 

jf  f^(x,0)^(x)dx  =  -  J^^^U{x,t)jlf{x,t)  . 
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Sfanilariy  for  Um  Auictioii  . «« baw 


X  !'»(*.«)♦(*)* = "X.,  'L,  • 

Thus,  by  (8.d),  w«  get 

Um  [  Uki*,0)^i*)dx  s  [  U(x,0)4(x)dx  V^(*)€i*(fi)  , 

4-^0  Jq  JQ 

Which  mesas  that  the  subsequence  1/4 (*.0)  — - -  Uix,0)  weakly  in  I»(n)  ,  and  so  we  have 

l^(x,0)  as  0  o.e.  on  fl  .  Similarly  we  can  show  that  o(*,0)  =  0  ,  o(x,0)  =  0  o.e.  on  0  ,  and 
all  the  initial  oonditkms  are  satisfied. 

Step  2.  lb  show  the  system  (8.5)  holds.  Indeed,  for  any  piecewise  bilinear  function  (fts), 

we  have 


V  ^4  €  <2i  (fls) 


(8.10) 


o.e.  t  €  /  .  So  for  any  ^(x,<)  €  CS*(ft)  x  6^(1)  ,  we  can  find  a  sequence  of  simple  abstract 
fonctions  ^4(*.0€  L^(I,Qtfi(^k)nHlia))  such  that  ^fc(x,0  converges  to  ^(x,t)  strongly 
in  if* (fix  7),  i.e. 

lim  f  (^4~^)*dxdl  =  0  (8*11) 

4-41  Jaxt 


Ita  /  (^-^)’<bdt  =  0 

4-4)  Jox/  OXi  axi 


tasl,2 


(8.12) 


For  any  » ,i,*  sa  1 ,2 ,  we  have 


(8.13) 


Prom  the  weak  *  convergence  of  the  sequence  04  ,  the  first  term  of  the  right  hand  side  of  (8.13)  will 
tend  to  zero  as  h  — »  0  .  For  the  second  term  of  the  right  hand  side  of  (8.13)  we  have  the  following 

inequality:  .  .  .  v  ,  11  o.  11 


1  Xx/**'  ^  1  -  I  ^  ^ 


lL»(nx/) 


<  ll^4vj|l,*(nx/)  1^4  ^ll««(Qx/) 

From  energy  estimates ,  we  know  that  ||  04,,  Hl^coxI)  5®  uniformly  bounded.  Meanwhile ,  from  the  strong 
H  *  convergence  of  04  to  0 ,  we  can  see  that  the  second  term  of  the  right  hand  side  of  (8.13)  will  also 

tend  to  zero .  Obviously 

lim  /  (/i,04)<^®*  =  L  »  =  1,2  . 

4-4)  Jnxr 
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St^  3.  Tb  show  that  the  limit  function  (^U^c^a)  still  catisfies  (8.7).  Let  (r,t/)  €  if  ,  then  we 
can  claim  that  there  exists  a  sequence  ( ,i% )  €  Xk  nich  that 


stron^y  in  ( fl  x  /)  . 


Indeed,  actually  we  can  construct  such  a  sequent  as  follows :  consider  a  sequence  of  moUiiiers  Jt  defined 
on  X  R  (  See  [1]  for  details  ) .  Then  we  have 


AW  =  ■/.(«)  =  ^ 

and 

Jt(x)dx  St  1  . 

Let  ft  ,  /  be  the  two  sets  such  that  ft  CC  ft  ,  /  CC  /  and  di8t(n,dft)  >  1  ,  di8t(/,d/)  >  1  . 
Define  r  (x)  s  v(x)  s  0  for  all  *  €  (R*  x  R)\(n  x  I).  For  any  s  €  H  x  7  and  c  <  1 ,  take 

*  L,  „  =  L  ,J€{y)T{z-y)dy  , 

-/R*xR  JQxI 


L  ^My)‘'{x-y)dy  =  f  .Jt{y)v(z-y)dy  . 

7R»xR  70x7 

Then  by  Lemma 2.18 of  [1],  we  have  T,(s),ye(s)  €  C^{(txi)  when  c  is  small  enough ,  and  t*  — ►  t, 
Vf  — »  V  strongly  in  x  /)  as  c  — »  0 .  On  the  other  hand,  from  the  Jenssen’s  Inequality  ( See,  for 
example  [7]  or  [16]  ),  for  any  convex  function  G  (u)  and  any  integrable  function  u(z  ,t)  over  Six  I , 
we  have 

or  equivalently 


where  m{Slx  I)  is  the  measure  of  the  set  0x7  and  C  is  an  arbitrary  constant.  So  for  our  gauge 
function  F(<T,a)  and  functions  Tt{z)  and  !/< (z)  ,  we  have 

7'(r*(z),*/.(z))  = /’(  /  Jt{y)T{z-  y)dy ,  f  7t(y)^(z  -  y)dy ) 

70x7  Jiix! 

-n - hi — Ti  L  .^(^tTn(Slxi)Jt{y)T{z-y),C,m{Slxi)Jt(y)i/{z-y))dy  . 

Ctin{Six  7)  70x7 
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'* 


as  - -T — 1-  ,  we  would  h»ve  CttnCllx  ^  y€ftx/.So  from  the 

km[UXi) 

convexity  oi  the  gange  function  F{t,v)  and  the  assumption :  F(0 ,0)  s  0  ,  we  get 

*  A, r(i  -  y)  +  (1  -  A,)0 ,  A, *<<*-»)  + (1  -  ■^)0 ) 

<  A,/’(T(s-y),i/(s-y))  +  (l- A«)/’(0,0)  =  A«  J’(r(s- y),i/(s- y))  , 
where  0  <  A,  ss  Cc  m  x  /)  J((y)  <  1 .  Hence  we  get 

^ (»«(*),*'€(*))<  L  Jt{M)F{r(x-y),v(z-y))dy<Ziil  Jciy)dy=Zo  , 

jQxl  i/R»xR 

which  implies  ( r*  ,»>« )  €  K  .  Because  ( r((y)  ,i/«(y)  )  €  C^{il  x  I)  ,  for  any  f  €  /  we  can  interp<date 
( ^c(y)  t  t'c(y) )  the  Gaussian  pdnts  by  the  basis  function  (^)  on  7\  and  construct  a  sequence  of 
simple  abstract  functions  (Tk,Vk)  *ach  that 

(TkiVk)  €  Kk  'i  ie  I 


and  {Tk,vk) — ►  strongly  in  L^{ClxI)  . 

Ihking  (t,i/)  *  (T)k,i/k)  in  (8.2),  we  get 


f  {ik-Ckk,Tk-Ok)dxdt -- {akyVk-ak)dxdt  <  0  , 

JQkJ 


which  can  be  written  in  the  form 


/  ((<fc-Cd*,Tfc)-(dA,yfc)ldxdt  -  /  (ak,ik)dxdt 
JOxl  JQxl 

The  weak*  convergence  of  the  sequence  (UkyCktCik)  uid  the  strong  convergence  of  the  sequence 
i'fkyVk)  will  ensure  that  the  first  term  tends  as  h  — »  0  to 


/  [{€-Cc)-{a,u)]dxdt  ; 

Jflxl 


the  second  term  is  equal  to 


f  {f,Uk)dxdt,  and  hence  tends  to  /  {/  ,U)dxdt 

Jflxl  Jflxl 

Because  the  initial  conditions:  <7’(z,0)  =  0,  a(z,0)  =  0,  the  third  term  is  equal  to 

{T)  Ili2(n)  I  » 
and 

inf  I  VCc,  (T)  >  |  VC  .(T)  , 

fclifSo  llL*(n)  ^  II  ®  {F)  11^(0)  • 
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0  with 


So  wo  ond  vp  M  h  —* 


lji-Ci.r)-ii,,)]-lj/,6)^.\[\VCam  +  II « (D  ||i.,„,]  <  0  . 


FVom  (8.5)  and  the  fact  that 


5ll«(r)|lli,o)  =  , 


WO  can  hnaOy  get  (8.7). 


Step  4.  To  show  that  (0,0)  €  A",  i.e.  Fio,Q)<Zo  o.e.  on  Qx/.Let  Pf^  :  R®xR"*i — ► 

X  It"*  be  the  projection  operator  on  the  convex  set 

a:  =  {(t,!/)  €  R®xR"  I  /’(t,i/)<Zo}  . 

For  the  convex  set  ,  we  can  ose  the  projection  theorem  and  get  that  for  any  (<r ,  a)  €  R^  x  R”*  and 

«®y  (r,i^)€  A” , 

Therefore  we  get 

For  any  (  €  /  ,  let  and  &h  be  the  piecewise  constant  functions  on  fl  such  that  on  each  rectangular 
dement  Tn  thdr  values  are  defined  by : 


(*fc)  »  ttk  -TT  a*  (**)  , 


(8.16) 


where  5n  is  the  area  of  the  rectangular  element  Tn  .  Then  from  the  fact  that  the  gauge  function 
F  ( o ,  a )  is  convex  and 


we  have 


Lf  =  1  . 

F  ( dfc ,  d* )  <  i;  ^  F  ( (Tfc(xfc) ,  Ofc(i*) )  <  E  ^  ^0  =  Zo  , 
/b=l  ibsl 


(dfc.dfc)  €  K 


o.e.  on  Q  X  / 


Consequently  we  have 


^  -  Pi(  °  t  ~  1  ^ 


(8.17) 
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fbr  tlM  function 
for  nny  fixed  1  € 


a 

a 


a 

a 


,  we  can  find  a  sequence  of  abstract  functibiis 


they  ate  piecewise  constant ,  and  as  h  — *  0 : 


Th 

0K 


such  that 


Tk 
0k  ^ 


L*lQxl) 


Now  we  can  rewrite  (8.17)  and  get 


Integrating  the  first  term  of  (8.19)  over  Q  x  7  , 


(8.18) 


(8.19) 


The  second  term  of  the  right  hand  side  of  the  abo^m  ec'^at  on  is  equal  to  zero,  because  ( rj^  ,0k )  are 
piecewise  constant  and  {ak,&k)  nre  defined  by  (8-16) .  So  from  (8.19)  we  get 


fVom  (8.18)  we  have ,  as  h  — »  0 , 


ff 

a 


a 

a 


-Pk 


t.e. 


a 

a 


-Pfi 


a.e.  on  fl  x  / 
and  so  (0,0)  €  K  ,ot  F(a,a)  <  Zq  a.e.  on  tlx  I . 


a 

a 


11 


<  0  , 


Step  5.  The  proof  of  the  uniqueness  follows  exactly  as  the  proof  of  the  Theorem  7.1  and  this  completes 
the  proof  of  the  theorem . 


For  the  solution  of  (8..5)~(8.8),  we  have  the  following  result : 

Proposition  8.1  The  solution  of  the  equation  system  (8.5)'^  (^•^)  solution  of  the  system 

(S.2)  '^(S.IO). 
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Ptoor:  Here  we  only  need  to  show  that  (8.2)  ~  (8.3)  imply  the  existence  of  the  scalar  function 

^(*>0  ^  0  such  that 


where 


{ 

{ 


€- C&  =  0 
a  =  0 

i-Ca  ^  ^  =  Xd^F 
a  =  -Xd^F 


a.e.  on  £ 


a.e.  on  V 


^s|(x,t)enx/  I  f’(<r,a)(*,<)  <  2o  or  ,0  <  0  } 

Ps|(x,<)€flx/  I  f’(o,ot)(x,t)  =  2b  and  ft,/’^d(x,t)  >  0 } 


(8.20) 


(8.21) 


Let  j  >  0  and  denote  by  (  xq  ,  to )  &  cube  of  side  length  S  centered  at  ( xo ,  to ) .  Then  for  any 
(  T ,  j/  )  €  fi’  ,  we  choose  (r,P)  such  that 

I  (f,i>)  =  (<T,a)  in  (ft  x /)\Q<  ( xo ,  to ) 

1  (^.^)  =  (’’♦*')  Q<(*o,to) 

Obviously  we  have  (f,^)  €  K  ,  and  (8.3)  yields 

f  [(€  -  C<T,T  -  <r)  +  (-d,i/ ~  o)]  dxdt  <  0  V(t,i/)€A'  • 

Using  the  Lebcsgue  OiiTerentiation  Theorem  (  See,  for  example,  [17] ),  we  get  as  6  — ►  0 

(€  -  Cd,r  -  o)  +  (-0,1/ -  o)  <  0  o.e.  on  fix/  V(t,i/)€A'  . 

Then  using  the  same  approach  as  we  used  in  the  step  3  of  the  proof  of  the  Theorem  7.1,  we  can  show  that 
there  exists  a  A(x,t)  >  0  such  that  (8.20)  and  (8.21)  hold .  □ 


Remark  8.1  Wc  can  easily  sec  that  this  method  can  be  extended  to  the  cases  where  the  displacement 
and  stress  functions  arc  approximated  by  piecewise  higher  order  polynomial  functions.  For  instance ,  if 
displacement  function  is  approximated  by  piecewise  or  functions ,  we  can  approximate  the 
corresponding  stress  function  by  piecewise  Qi  functions ,  and  in  this  case ,  nine  Gaussian  points  will  be 
used  for  the  constitutive  equations . 


Remark  8.2  Unfortunately  this  semi-discrete  approximation  method  can  only  be  directly  appUed  to 
the  cases  where  the  partition  of  the  domain  fl  only  contains  rectangular  and  parallelogram  elements . 
However ,  as  we  will  see  in  the  foUowing  section ,  for  cases  where  the  partition  of  the  domain  contains 
some  other  type  of  elements ,  similar  approach  can  be  used  to  develop  3ome  higher  order  methods . 
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9  Generalization  to  the  Triangular  Elements 


It  it  detirable  to  find  a  method  which  can  be  applied  to  triangular  dements  and  that  alto  hat  the 
tame  properties  at  rectangular  or  paralldogram  dements .  Then  we  would  be  able  to  treat  cases  where 
SI  is  any  pdygonal  domain .  In  this  section ,  we  will  briefly  discuss  a  method  of  using  piecewise  higher 
order  polynomials  over  triangular  dements  . 


Consider  the  natural  triangle  of  area  A  as 
shown  in  Figure  9.1 ,  where  the  natural  coordi¬ 
nates  (0,^,7)  are: 


where  Ai,  A^,  M  and  A  are  the  area  of  each 
triangles .  Hence  we  have : 

0  =  a-i-i8-i-7-l  . 

Now  consider  another  standard  triangle  on 
the  X  -  Y  plane  as  shown  in  Figure  9.2 ,  and 
take  any  point  ( z ,  y )  inside  the  triangle . 


a 


Figure  9.1  Natural  Triangle 
and  Coordinates 


Figure  9.2  Standard  Triangle 
on  the  X-Y  Plane 


Then  a  simple  calculation  shows : 

A-^/Z 

Ai=y 

i42  =  -  “('N/Sz  +  y  — >/3) 

if 

A3  =  i(V3z-y-|-v/3) 

Therefore  we  get  a  one-to-one  correspon¬ 
dence  between  the  X-Y  coordinate  system  and 
the  natural  coordinate  system  given  by  the  fol¬ 
lowing  mapping  functions : 


7  =  + 


(9.1) 
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Therefore: 


with  dA  s  2Adad0  . 


I  fi*,y)dxdy=  I  f(oi,fi,-i)dA 


Some  symmetrical  Gaussian  quadrature  rules  were  discussed  in  [5] ,  [6]  and  [9]  ,  where  an  integration 
is  performed  by  a  Gaussian  quadrature  rule  of  the  form : 

/(Q,)9,7)di4  =  i4  .  (9.2) 

ml 

In  (9.2) ,  (oti ,  ,7i)  are  the  natural  coordinates  of  the  tth  Gaussian  point,  a;,-  the  corresponding 

Gaussian  weight  and  ng  is  the  number  of  Gaussian  points  used  in  the  rule  . 


Table  9.1  Symmetrical  Gaussian  Quadrature  Rules 


p  value 

»g 

weight 

alpha 

beta 

gamma 

1 

1 

1.0 

1/3 

1/3 

1/3 

2 

3 

1/3 

2/3 

1/6 

1/6 

im 

1 

0.223381590 

0.108103018 

0.445948491 

0.445948491 

0.109951744 

0.816847573 

0.091576214 

0.091576214 

5 

1 

0.225000000 

0.333333333 

0.333333333 

0.333333333 

0.132394152 

0.059715872 

0.470142064 

0.470142064 

0.125939181 

0.797426985 

0.101286507 

0.101286507 

6 

12 

0.116786276 

0.501426509 

0.249286745 

0.249286745 

0.050844906 

0.873821971 

0.063089014 

0.063089014 

0.082851076 

0.053145050 

0.310352451 

0.636502499 

8 

16 

0.144315677 

0.333333333 

0.333333333 

0.333333333 

0.095091634 

0.081414823 

0.459292588 

0.459292588 

0.103217371 

0.658861384 

0.170569308 

0.170569308 

0.032458498 

0.898905543 

0.050547228 

0.050547228 

0.027230314 

0.008394777 

0.263112830 

0.728492393 

9 

19 

0.097135796 

0.333333333 

0.333333333 

0.333333333 

0.031334700 

0.020634962 

0.489682520 

0.489682520 

0.077827541 

0.125820817 

0.437089591 

0.437089591 

0.079647739 

0.623592929 

0.188203536 

0.188203536 

0.025577676 

0.910540973 

0.044729513 

0.044729513 

0.043283539 

0.036838412  |  0.221962989 

0.741198599 

A  table  of  Gaussian  quadrature  rules  of  the  form  (9.2)  for  the  polynomials  of  degree  1-20  was 
Usted  in  [6] .  Different  from  the  Gaussian  quadrature  rules  over  rectangles ,  the  Gaussian  quadrature 
rules  presented  in  [6]  have  some  negative  weight  coefficients  u;,-  or  some  Gaussian  points  not  lying 
inside  the  triangle  A  for  some  p  values .  In  our  case  we  cannot  use  Gaussian  quadrature  rules  with 
negative  weights ,  because ,  as  we  have  seen  m  the  previous  discussion ,  the  positiveness  of  the  weight 
coefficients  plays  a  very  important  role  in  the  proofs  of  many  theorems.  We  also  cannot  use  those 
Gaussian  quadrature  rules  with  some  Gaussian  points  lying  outside  the  triangle ,  because  we  have  to  use 
extrapolation  to  determine  the  stress  and  hardening  parameter  functions.  So  in  Table 9.1  we  only  list 
the  rules  with  positive  weights  and  all  Gaussian  points  inside  the  triangle  A  . 
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In  Ti&ble  9.1,  for  each  weight  if  we  have  a  =  =  7  ,  then  there  is  only  one  Gaussian  point 
corresponding  to  that  wdght  coefficient .  If  we  have  a  /  ^  =  7  ,  then  there  are  C3  =  3  Gaussian 
pdnts  corresponding  to  that  weight  coefficient  .  Finally  if  we  have  a  ^  ^  7  ,  then  there  are  C3  s  6 

Ganssian  points  corresponding  to  that  weight  coefficient .  Here  we  only  list  the  weights  and  the  location 
of  the  Gaussian  points  for  the  p  values  up  to  9 ,  for  more  data  with  9  <  p  <  20  see  [6] . 

Now  we  can  discuss  the  selection  of  the  finite  dimensional  spaces  for  the  displacement  function 
1/  (x ,  t)  and  the  stress  and  hardening  parameter  functions  o  (x ,  t ) ,  a  (x ,  t) .  The  strata  is  as  follows : 
the  finite  dimensional  spaces  for  U  ^  a  and  a  should  be  selected  in  such  a  way  that 

i )  .  The  stress  and  hardening  parameter  functions  a ,  a  will  be  uniquely  determined  by  their  values 
at  the  ng  Gaussian  points .  (  Assume  here  we  want  to  use  ng  Gaussian  Points ) 

ii) .  The  degree  of  the  polynomials  of  functions :  and  a  •  e(U)  must  be  less  than  or  equal  to 

the  p  value  listed  in  Table  9.1  corresponding  to  the  number  ng  . 

Table  9.2  gives  a  list  of  a  possible  selections  of  the  spaces  for  U ,  a  and  a  .  Here  we  always  select 
the  same  spaces  for  a  and  a  . 

Table  9.2  A  Possible  Selection  of  the  Spaces  for  U ,  a  and  a  . 


ng 

D.O.P. 

Space  for  U 

Space  for  0  and  a 

1 

1 

Pi 

Po 

3 

2 

Pj 

Pi 

6 

4 

Ps 

Pi 

For  the  complete  polynomial  subspaces ,  we  can  only  get  these  three  possible  combinations .  How¬ 
ever,  it  is  not  necessary  to  always  use  complete  polynomial  subspaces  for  the  displacement  ,  stress 
and  hardening  functions.  As  a  matter  in  fact,  by  using  a  Gaussian  quadrature  rule  based  on  some 
non-complete  polynomial  subspaces ,  we  can  also  get  some  higher  order  spaces  for  U ,  a  and  a  . 

It  is  also  not  necessary  to  always  choose  piecewise  polynomial  functions  as  the  subspace  for  the 
stress  and  hardening  parameter  functions .  Since  we  don't  require  that  a ,  a  are  functions  over  the 
domain  fl  .  Actually  we  can  select  any  n  dimensional  subspace  Sn  €  over  each  rectangle  or 

triangle  element ,  if  this  subspace  has  the  following  properties : 

i )  .  The  stress  and  hardening  parameter  functions  in  this  subspace  will  be  uniquely  determined  by  their 
values  at  some  sampling  points  Xk  k  =  \,...,n. 

ii )  .  There  is  a  quadrature  rule  over  the  rectangle  or  triangle  element  with  non-negative  weights  w* 
such  that 

/» 

f{x,y)dxdy=  X^u;fc/(x*,y*) 
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for  all  the  /  (x ,  y )  of  the  forms : 

(e(f/(x,ir)))’,  €(U(x,y)).g(x,y),  {g{x,y)f 

where  U  {x  ,y)  is  any  displacement  function  in  its  finite  dimensional  space  and  g{x  ,y)  is  any  function 
in  the  space  Sn  • 

Finally  we  would  lilce  to  indicate  that  the  ideal  of  using  Gaussian  point  may  also  be  applied  to  the 
cases  where  the  partition  of  the  domain  contains  some  curved  elements .  The  convergence  result  we  get 
here  is  only  in  the  weak  L^(il  x  I)  sense,  which  is  obviously  not  good  enough  for  real  computational 
purpose .  However ,  we  can  prove  that  for  some  special  constitutive  laws ,  such  as  bilinear  isotropic  or 
bilinear  kinematic  hardening  laws ,  we  can  actually  get  strong  L^(ClxI)  convergence .  For  more  general 
constitutive  laws ,  the  strong  convergence  is  still  achievable ,  if  we  use  some  extra  assumptions  on  the  the 
body  traction  functions . 
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